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Abstract. We report a measurement of the time-dependent ratio of D" — > K*n^ to £)" — » K^n* decay rates 
in Z)*^ -tagged events using 1.0 fb"' of integrated luminosity recorded by the LHCb experiment. We mea- 
sure the mixing parameters x'^ = (-0.9 ± 1.3) x 10"'', y' = (7.2 ± 2.4) x 10"^ and the ratio of doubly-Cabibbo- 
suppressed to Cabibbo-favored decay rates Ro = (3.52 ± 0.15) x 10"^. The result excludes the no-mixing hy- 
pothesis with a probability corresponding to 9.1 standard deviations and represents the first observation of 
charm mixing from a single measurement. 



1 Introduction 

Quantum-mechanical mixing between neutral meson par- 
ticle and antiparticle flavour eigenstates provides impor- 
tant information about electroweak interactions and the 
Cabibbo-Kobayashi-Maskawa matrix, as well as the vir- 
tual particles that are exchanged in the mixing process 
itself. For this reason, the mixing of neutral mesons is 
generally considered a powerful probe to discover physics 
beyond the standard model. Meson-antimeson mixing has 
been observed in the K°-K^^ \ 1 1, [2] and B'j-B? S 

systems, all with rates in agreement with standard model 
expectations. Evidence of mixing in the charm system 
has been reported by three experiments using different 

decay channels |4-9|. Only the combination of these 
measurements provides confirmation of Z)" - D" mixing 
with more than 5cr significance ifTOl . While it is accepted 
that charm mixing occurs, a clear observation of the phe- 
nomenon from a single measurement is needed to establish 
it conclusively. 

Thanks to the large charm production cross-section 
available in pp collisions at = 7 TeV and to its flexible 
trigger on hadronic final states, the LHCb experiment IfTTI 
collected the world's largest sample of fully reconstructed 
hadronic charm decays during the 2011 LHC run. In the 
following we present the first search for - mixing 
using this data sample, which corresponds to l.Ofb ' of 
integrated luminosity. More details on the analysis can be 
found in Ref. I.12J . 

2 IVIeasurement of charm mixing in thie 

^ K^n chiannel 

In this analysis a search for - mixing is performed 
by measuring the time-dependent ratio of D" — > K^n^ 
to D" ^ K decay ratesj^ The D" flavour at pro- 
duction time is determined using the charge of the soft 

' Charge-conjugated modes are implied unless otherwise stated. 



(low-momentum) pion, n'^, in the strong D*^ D^^t 
decay. The D*^ — » D"(— > K^7:^)7r^ process is referred 
to as right-sign (RS), whereas the D*+ D^{^ K^n')n^ 
is designated as wrong-sign (WS). The RS process is 
dominated by a Cabibbo-favored (CF) decay amplitude, 
whereas the WS amplitude includes contributions from 
both the doubly-Cabibbo-suppressed (DCS) Z)° -> K^n~ 
decay, as well as D" - D" mixing followed by the fa- 
vored D" — > K^n^ decay. In the limit of small mixing 
(|x|,|i/| 1), and assuming negligible CP violation, the 
time-dependent ratio, R(t), of WS to RS decay rates is ap- 
proximated by 

R{t) ^ Rd + yfR^ y' - + — ^("J' (1) 

where f/r is the decay time expressed in units of the aver- 
age D'^ lifetime t, R^is the ratio of DCS to CF decay rates, 
x' and y' are the mixing parameters rotated by the strong 
phase difference between the DCS and CF amplitudes. In 
case of no mixing, x' - y' - 0, the WS/RS ratio would be 
independent of decay time. 

2.1 Outline of the analysis 

The candidate reconstruction exploits the capabilities of 
the silicon vertex locator (VELO), to identify the dis- 
placed two-track vertices of the decay products with 
decay-time resolution Af a: O.lr; of the tracking system, 
which measures charged particles with momentum reso- 
lution A/?//? that varies from 0.4% at 5 GeV/c to 0.6% at 
100 GeV/c, corresponding to a typical mass resolution of 
approximately 8 MeV/c^ for a two-body charm-meson de- 
cay; and of the ring imaging Cherenkov detectors, which 
are used to distinguish between pions and kaons and to 
suppress the contamination from misidentified two-body 
charm decays in the sample. 

We reconstiTict approximately 3.6 x 10"* WS and 8.4 x 
10^ RS decays as shown in figure 111 where the M(D"7r^) 
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Figure 1. Time-integrated D^n^ mass distributions for the selected RS ^ K n* (left) and WS D" K*n (right) candidates with 
fit projections overlaid. The bottom plots show the normalized residuals between the data points and the fits. 



distribution for the selected RS and WS candidates is fitted 
to separate the D*^ signal component, with a mass reso- 
lution of about 0.3 y[€^lc^, from the background compo- 
nent, which is dominated by associations of real D*' decays 
and random pions. Similar fits are used to determine the 
signal yields for the RS and WS samples in thirteen D" 
decay time bins, chosen to have a similar number of can- 
didates in each bin. The shape parameters and the yields of 
the two components, signal and random pion background, 
are left free to vary in the different decay time bins. We 
further assume that the M{D^jil) signal shape for RS and 
WS decay are the same and therefore first perform a fit 
to the more abundant and cleaner RS sample to determine 
the signal shape and yield, and then use those shape pa- 
rameters with fixed values when fitting for the WS signal 
yield. The signal yields from the thirteen bins are used to 
calculate the WS/RS ratios and the mixing parameters are 
determined in a binned fit to the observed decay-time 
dependence. 

2.2 Systematic uncertainties 

Since WS and RS events are expected to have the same 
decay-time acceptance and MiD^n^) distributions, most 
systematic uncertainties affecting the determination of the 
signal yields as a function of decay time cancel in the 
ratio between WS and RS events. Residual biases from 
non-cancelling instrumental and production effects, such 
as asymmetries in detection efficiencies or in production, 
are found to modify the WS/RS ratio only by a relative 
fraction of 0(10"'*) and are neglected. Uncertainties in the 



distance between VELO sensors can lead to a bias of the 
decay-time scale. The effect has been estimated to be less 
than 0.1% of the measured time and translates into relative 
systematic biases of 0.1% and 0.2% on y' and x'^, respec- 
tively. At the current level of statistical precision, such 
small effects are negligible. 

The main sources of systematic uncertainty are those 
which could alter the observed decay-time dependence of 
the WS/RS ratio. Two such sources have been identi- 
fied: (1) D mesons from fo-hadron decays, and (2) peak- 
ing backgrounds from charm decays reconstructed with 
the wrong particle identification assignments. These ef- 
fects, discussed below, are expected to depend on the true 
value of the mixing parameters and are accounted for in 
the time-dependent fit. 

A contamination of charm mesons produced in b- 
hadron decays (secondary D decays) could bias the time- 
dependent measurement, as the reconstructed decay time 
is calculated with respect to primary vertex, which, for 
these candidates, does not coincide with the pro- 
duction vertex. When the secondary component is not 
subtracted, the measured WS/RS ratio can be written 
as R{t)[\ - AB(f)], where R{t) is the ratio of promptly- 
produced candidates according to Eq. ([T]i, and Ab(0 is 
a time-dependent bias due to the secondary contamina- 
tion. Since R{f) is measured to be monotonically non- 
decreasing 1 10 1 and the reconstructed decay time for sec- 
ondary decays overestimates the true decay time of the D" 
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Figure 2. Background-subtracted (IP) distributions for RS £)" candidates in two different decay-time bins. The dashed line indicates 
the cut used in the analysis; the solid histograms represent the estimated secondary components. 
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Figure 3. Measured fraction of secondary decays entering the fi- 
nal RS sample as a function of decay time (points), with overlaid 
the projection of a fit to a sigmoid-like function (solid line). 



meson, it is possible to bound Agd), for all decay times, as 



< AbH) < /f (0 



1 - 



(2) 



where /^^(f) is the fraction of secondary decays in the 
RS sample at decay time t. In this analysis most of the 
secondary D decays are removed by requiring the ;t'^(IP) 
of the D° to be smaller than 90 A residual (2.7 ± 0.2)% 



contamination survives. To include the corresponding sys- 
tematic uncertainty, we modify the fitting function for the 
mixing hypothesis assuming the largest possible bias from 
equation (|2]i. The value of f^^(t) is constrained, within 
uncertainties, to the measured value, obtained by fitting 
the ;^^(IP) distribution of the RS D'^ candidates in bins 
of decay time (see figure [2]). In this fit, the promptly- 
produced component is described by a time-independent 
X^(iP) shape, which is derived from data using the candi- 
dates with t < 0.8t. The ;\f^(IP) shape of the secondary 
component (represented by the solid histograms in fig- 
ure |2]i, and its dependence on decay time, is also deter- 
mined from data by studying the sub-sample of candidates 
that are reconstructed, in combination with other tracks in 
the events, as B ^ £>*(3);r, B D*nX or B ^ D°ijX. 
The measured value of fg^(t) is shown in figure |3] As- 
suming the maximum bias could induce an over-correction 
which results in a shift in the estimated mixing parame- 
ters. We checked on pseudo-experiments, before fitting 
the data, and then also on data that such a shift is always 
much smaller than the corresponding increase in the un- 
certainty when the secondary bias is included in the fit. 

Peaking background in M(D°7r^), that is not accounted 
for in our mass fit, arises from D*^ decays for which the 
correct soft pion is found but the is partially recon- 
structed or misidentified. This background is suppressed 
by the use of tight particle identification and two-body 
mass requirements. From studies of the events in the D" 
mass sidebands, we find that the dominant source of peak- 
ing background leaking into our signal region is from RS 
events which are doubly misidentified as a WS candidate; 
they are estimated to constitute (0.4 ± 0.2)% of the WS 
signal. From the same events, we also derive a bound 



■^Thex^i^P) is defined as the difference between the^^ of the primary 
vertex reconstructed with and without the considered particle, and is a 



measure of consistency with the hypothesis that the particle originates 
from the primary vertex. 



xlO" 



Q 0.25 : 



C/3 



o 



0.2 



0.15 



0.05 







I— p-i— I— I— p-i— I— I— p-i— I— I— p-i— I— 
Data LHCb 

Linear fit with 
68% C.L. band 

Constant fit 




4 



20 



^ 7 
6.5 

6 
5.5 

5 
4.5 

4 
3.5 

3 



xlO" 







_ 1 1 1 1 1 1 1 1 1 1 1 

: • Data 


1 1 1 1 1 1 1 ^ 


r — Mixing fit 




r No-mixing fit 




: I 

: 




LHCb J 

~ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 / /i 1 1 " 



20 
tlx 



Figure 4. Decay-time evolution of the number of doubly 
misidentified RS events observed in the £)" mass sidebands of the 
WS sample normalized to the RS signal yield. The solid (dashed) 
line is the result of a fit assuming linear (constant) decay-time de- 
pendence. 



Figure 5. Decay-time evolution of the ratio, R, of WS 
£)" -» K^tC to RS £)" -> K'n* yields (points) with the projec- 
tion of the mixing allowed (solid line) and no-mixing (dashed 
line) fits overlaid. 



on the possible time dependence of this background (see 
figure |4|i, which is included in the fit in a similar manner 
to the secondary background. Contamination from peak- 
ing background due to partially reconstructed D° decays is 
found to be much smaller than 0. 1 % of the WS signal and 
neglected in the fit. 

2.3 Results 

Figure |5] shows the observed decay -time evolution of the 
WS to RS ratio, with the projection of the best fit result 
overlaid (solid line). The estimated values of the parame- 
ters Rd, y' and x are listed in table [l] 

Table 1. Results of the time-dependent fit to the data. The 
uncertainties include statistical and systematic sources; ndf 
indicates the number of degrees of freedom. 



Fit type Or^/ndf) 


Parameter 


Fitresuh(10"^) 


Mixing (9.5/10) 


Rd 


3.52 ±0.15 




1/ 


7.2 + 2.4 






-0.09 + 0.13 


No mixing (98.1/12) 


Rd 


4.25 + 0.04 
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Figure 6. Estimated confidence-level (CL) regions in the {x'^,y') 
plane for 1 - CL = 0.317 (Icr), 2.7 x 10"^ (3cr) and 5.73 x 10"' 
(5cr). Systematic uncertainties are included. The cross indicates 
the no-mixing point. 



To evaluate the significance of this mixing result we 
determine the change in the fit when the data are de- 
scribed under the assumption of the no-mixing hypothesis 
(dashed line in figure [5]). Under the assumption that the 

difference, Ay^, follows a x^ distribution for two de- 
grees of freedom, Ay^ = 88.6 corresponds to a p-value of 
5.7 X 10"^°, which excludes the no-mixing hypothesis at 
9.1 standard deviations. This is also illustrated in figure |6] 



where the Icr, 3cr and 5cr confidence regions for x'^ and 
are shown. 



3 Conclusions 

We measure the decay time dependence of the ratio be- 
tween — > K^n'' and D" — > K^n* decays using 1.0 fb"' 



of data and exclude the no-mixing hypothesis at 9.1 stan- 
dard deviations. This is the first observation of D° -Z)° os- 
cillations in a single measurement. The measured values 
of the mixing parameters are compatible with and have 
better precision than those from previous measurements 
H mini, as shown in table |2] 

Table 2. Comparison of our result with recent measurements 
from other experiments. The uncertainties include statistical and 
systematic components. 



Exp. 


Rd (10-^) 


y' i^o-') 


(10-^) 


LHCb 2011 


3.52 + 0.15 


7.2 ± 2.4 


-0.09 ±0.13 


BaBarim 


3.03 + 0.19 


9.7 ± 5.4 


-0.22 ± 0.37 


Belle ifBl 
CDFUa 


3.64 + 0.17 
3.04 ± 0.55 


8.5 ±7.6 


-0.12 ±0.35 



The result is still statistically limited and can be im- 
proved thanks to the already available 2.1 fb ' sample of 
pp collisions recorded by LHCb at -v/i = 8 TeV during 
2012. This additional sample will also allow to search 
for CP violation in charm mixing with sensitivities never 
achieved before. 
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